Similarity- Transformed Chiral NN+3N Interactions for the Ab Initio Description of 12 C and ie O 
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We present first ab initio no-core shell model (NCSM) calculations using similarity renormalization group 
(SRG) transformed chiral two-nucleon (NN) plus three-nucleon (3N) interactions for nuclei throughout the p- 
shell, particularly 12 C and 16 0. By introducing an adaptive importance truncation for the NCSM model space 
and an efficient yr-coupling scheme for the 3N matrix elements, we are able to surpass previous NCSM studies 
including 3N interactions regarding, both, particle number and model-space size. We present ground and excited 
states in 12 C and 16 for model spaces up to iV max =12 including full 3N interactions, which is sufficient to obtain 
converged results for soft SRG-transformed interactions. We analyze the contributions of induced and initial 
3N interactions and probe induced 4N terms through the sensitivity of the energies on the SRG flow parameter. 
Unlike for light p-shell nuclei, SRG-induced 4N contributions originating from the long-range two-pion terms 
of the chiral 3N interaction are sizable in 12 C and 16 0. 
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Nuclear Hamiltonians constructed within chiral effective 
field theory (EFT) provide a systematic link between nuclear 
structure physics and low-energy quantum chromodynamics 
(QCD). It is a supreme challenge for modern nuclear theory 
to exploit this link and to apply these interactions consistently 
in ab initio nuclear structure calculations for a wide range of 
nuclei and observables. This is vital to provide robust QCD- 
based predictions, e.g. for light exotic nuclei, to constrain ap- 
proximate nuclear structure approaches, and to understand the 
relevant QCD mechanisms driving nuclear structure phenom- 
ena. At present the most advanced calculations beyond the 
few-body domain use chiral two-nucleon (NN) interactions at 
N 3 LO lUl] and three-nucleon (3N) interactions at N 2 LO 0] in 
the ab initio no-core shell model (NCSM) ft^. A first mile- 
stone was the study of the spectroscopy of mid-p-shell nu- 
clei in moderate model spaces using a Lee-Suzuki (LS) trans- 
formed Hamiltonian 0], which proved the predictive power 
of chiral Hamiltonians and the need to include the 3N interac- 
tion consistently. Recently, light p-shell nuclei were studied in 
the NCSM employing a consistent similarity renormalization 
group(SRG) transformation of the chiral NN+3N Hamilto- 
nian JH101. The SRG transformation provides a model-space 
independent Hamiltonian with superior convergence proper- 
ties that can be used universally in a variety of many-body 
approaches 

In this Letter we present the first ab initio calculations of 
nuclei throughout the whole p-shell including 12 C and I6 us- 
ing SRG-transformed chiral NN+3N interactions. Through a 
combination of conceptual and computational developments 
we are able to extend the range of previous NCSM studies us- 
ing full 3N interactions to significantly larger model spaces 
and particle numbers. 

SRG-Transformed NN+3N Interactions. A crucial step 
for NCSM calculations beyond the lightest isotopes is the uni- 
tary transformation of the initial Hamiltonian in order to im- 
prove the convergence behavior with respect to the size of 
the many-body model space. In addition to the LS similar- 



ity transformation, which is tailored to decouple the NCSM 
model space from the excluded space, several model-space 
independent unitary transformations, e.g. the unitary correla- 
tion operator method (UCOM) and the SRG, have been intro- 
duced m. Here, we focus on the SRG, mainly because of 
its simplicity and flexibility. 

In the SRG framework the unitary transformation of an op- 
erator, e.g. the Hamiltonian, is formulated in terms of a flow 
equation -j^H a = [rj a ,H a ] with a continuous flow parameter 
a. The initial condition for the solution of this flow equation 
is given by the 'bare' chiral Hamiltonian. The physics of the 
SRG evolution is governed by the anti-hermitian generator r\ a . 
A specific form widely used in nuclear physics tzl |8|] is given 
by r\ a = tn 2 N [T mt ,H a ], where is the nucleon mass and 
Tint = T — T cm is the intrinsic kinetic energy operator. This 
generator drives the Hamiltonian towards a diagonal form in 
a basis of eigenstates of the intrinsic kinetic energy. 

Along with the pre-diagonalization of the Hamiltonian, 
which is the reason for the transformation in the first place, 
the SRG induces many-body operators beyond the particle- 
rank of the initial Hamiltonian. Only if all the induced terms 
up to the A -body level are kept, the transformation is unitary 
and the spectrum of the Hamiltonian in an exact A-body cal- 
culation is unchanged and independent of the flow parameter 
a. In practice we have to truncate the evolution at a particle 
rank n < A, thus violating formal unitarity. In this situation we 
can use the flow-parameter a as a diagnostic tool to quantify 
the contribution of omitted beyond-«-body terms. 

Whereas the SRG transformation at two-body level has 
been used for some time the solution of the evo- 

lution equations at three-body level was demonstrated only 
recently J5l0]. In view of the application in the NCSM it is 
convenient to solve the flow equation for the three-body sys- 
tem using a harmonic-oscillator (HO) Jacobi-coordinate basis 
lUlll . The intermediate sums in the 3N Jacobi basis are trun- 
cated at Mnax = 40 for channels with J < 5/2 and ramp down 
linearly to Mn ax = 24 for J > 13/2. Based on this and the cor- 
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responding solution of the flow equation in two-body space 
(using either a partial-wave momentum or harmonic-oscillator 
representation) we extract the irreducible two- and three-body 
terms of the Hamiltonian for the use in A-body calculations. 

We have made major technical improvements regarding 
the SRG transformation, reducing the computational effort by 
three orders of magnitude compared to Ref. H, e.g., by us- 
ing a solver with adaptive step-size and optimized matrix op- 
erations. Furthermore, we have developed a transformation 
from 3N Jacobi matrix elements to a /r-coupled representa- 
tion with a highly efficient storage scheme, which allows us 
to handle 3N matrix-element sets of unprecedented size. A 
detailed discussion of these aspects is presented elsewhere. 

Importance-Truncated NCSM. Based on the SRG- 
evolved Hamiltonian we treat the many-body problem in the 
NCSM, i.e., we solve the large-scale eigenvalue problem 
of the Hamiltonian, represented in a many-body basis of 
HO Slater determinants truncated w.r.t. the maximum HO 
excitation energy N max nQ.. In order to cope with the factorial 
growth of the basis dimension with N max and particle number 
A, we use the importance-truncation (IT) scheme introduced 
in Refs. Il2lll3ll . The IT-NCSM uses an importance measure 
k v for the individual basis states |O v ) derived from many -body 
perturbation theory and retains only states with \k v \ above 
a threshold K mm in the model space. Through a variation 
of the threshold and an a posteriori extrapolation K m \ n — > 
the contribution of discarded states is recovered. We use 
the sequential update scheme discussed in Ref. [13], which 
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connects to the full NCSM model space and thus the exact 
NCSM results in the limit of vanishing threshold. In the 
following we always report threshold-extrapolated results 
including an estimate for the extrapolation uncertainties. For 
the present application we have extended the IT-NCSM to 
include full 3N interactions. Using the /r-coupled 3N matrix 
elements we are able to perform calculations up to N nYdx = 12 
or 14 for all p-shell nuclei with moderate computational 
resources. Due to the /r-coupling, we can keep all 3N matrix 
elements in memory using a fast on-the-fly decoupling. 

Ground-State Energies. We first focus on IT-NCSM cal- 
culations for the ground states of 4 He, 6 Li, 12 C, and I6 us- 
ing SRG-transformed chiral NN+3N interactions. Through- 
out this work we use the chiral NN interaction at N 3 LO of 
Entem and Machleidt Htl and the 3N interaction at N 2 LO J3 
with low energy constants determined from the triton binding 
energy and /?-decay half-live H 1 511 - In order to disentangle the 
effects of the initial and the SRG-induced 3N contributions, 
we consider three different Hamiltonians. (1) NN only: start- 
ing from the chiral NN interaction only the SRG-evolved NN 
contributions are kept. (2) NN+3N -induced: starting from the 
chiral NN interaction the SRG-evolved NN and the induced 
3N terms are kept. (3) NN+3N-full: starting from the chiral 
NN+3N interaction the SRG-evolved NN and all 3N terms 
are kept. For each Hamiltonian we assess the dependence of 
the observables, here the ground-state energies, on the flow- 
parameter a. We use the five values a = 0.04 fm 4 , 0.05 fm 4 , 
0.0625 fm 4 , 0.08 fm 4 , and 0.16 fm 4 , which correspond to mo- 
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FIG. 1 : (color online) IT-NCSM ground-state energies for 4 He and 
6 Li as function of N max for the three types of Hamiltonians (see col- 
umn headings) for a range of flow parameters: a = 0.04 fm 4 (•), 
0.05 fm 4 (♦), 0.0625 fm 4 (a), 0.08 fm 4 (■), and 0.16fm 4 (★). Error 
bars indicate the uncertainties of the threshold extrapolations. The 
bars at the right-hand-side of each panel indicate the results of expo- 
nential extrapolations of the individual _/V raax -sequences (see text). 



mentum scales A = a' 114 = 2.24 fm" 1 , 2.11 fm" 1 , 2fm _1 , 
1.88 fm" 1 , and 1.58 fm" 1 , respectively. For extrapolations to 
infinite model space, N max — » oo, we use simple exponen- 
tial fits based on the last 3 or 4 data points. The extrapolated 
energy is given by the average of the two extrapolations, the 
uncertainty by the difference. 

The ground-state energies obtained in IT-NCSM calcula- 
tions for 4 He and 6 Li with the three Hamiltonians are sum- 
marized in Fig. Q] Analogous calculations in the full NCSM 
for the same SRG-evolved initial Hamiltonian have been pre- 
sented in Ref. H for 4 He and in Ref. 01 for 6 Li. We have 
cross-checked our results with Refs. and found excel- 

lent agreement. 

The first and foremost effect of the SRG transformation 
is the acceleration of the convergence of NCSM calculations 
with 7V max . With increasing a the convergence is systemati- 
cally improved for all three versions of the Hamiltonian. With 
the initial Hamiltonian, i.e. a = 0, even the large model spaces 
we use here are not sufficient to even obtain meaningful ex- 
trapolations, with the exception of the tightly-bound 4 He. 

For the NN-only Hamiltonian Fig. Q] shows a clear a- 
dependence of the extrapolated ground-state energies for 4 He 
and 6 Li, hinting at sizable SRG-induced 3N contributions. 
When including those induced 3N terms, i.e. when using the 
NN+3N-induced Hamiltonian, the extrapolated ground-state 
energies are shifted significantly and become (^-independent 
within the uncertainties of the 7V max -extrapolation. Thus, in- 
duced contributions beyond the 3N level originating from the 
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FIG. 2: (color online) IT-NCSM ground-state energies for 12 C and 
le O as function of N mwi for the three types of Hamiltonians and a 
range of flow parameters (for details see Fig. Q]. 



initial NN interaction are negligible in the ff-range considered 
here, indicating that the NN+3N-induced Hamiltonian is uni- 
tarily equivalent to the initial NN Hamiltonian. The extrapo- 
lated ground-state energies for different a are summarized in 
TabJ 

By including the initial chiral 3N interaction, i.e., by using 
the NN+3N-full Hamiltonian, the ground-state energies are 
lowered and are in good agreement with experiment for both, 
4 He and 6 Li. As for the NN+3N-induced there is no sizable 
ff-dependence in the range considered here. We conclude that 
induced 3N terms originating from the initial NN interaction 
are important, but that induced 4N (and higher) terms are not 
relevant for light p-shell nuclei, since the ground-state ener- 
gies obtained with the NN+3N-induced and the NN+3N-full 
Hamiltonian are practically a-independent. 

This picture changes if we consider nuclei in the upper p- 
shell. In Fig. [2] we show the first accurate ab initio calcula- 



TABLE I: Summary of ,/V max -extrapolated IT-NCSM ground-state en- 
ergies in MeV for Tiil = 20 MeV (see text). 





a [fm 4 ] 


4 He 


"Li 


12 C 


16q 


NN 
only 


0.05 

0.0625 

0.08 


-28.08(2) 
-28.25(1) 
-28.38(1) 


-31.5(2) 
-31.8(1) 
-32.2(1) 


-99.1(6) 

-101.4(3) 

-103.7(2) 


-161.0(2) 
-164.9(6) 
-170.2(4) 


NN+ 
3N-ind. 


0.05 

0.0625 

0.08 


-25.33(1) 
-25.34(1) 
-25.34(1) 


-27.7(2) 
-27.6(2) 
-27.6(1) 


-76.9(2) 
-77.2(1) 
-77.4(2) 


-119.5(3) 
-119.7(6) 
-119.5(2) 


NN+ 
3N-full 


0.05 

0.0625 

0.08 


-28.45(3) 
-28.45(1) 
-28.46(1) 


-31.8(2) 
-31.8(1) 
-31.8(1) 


-96.1(4) 
-96.8(3) 
-97.6(1) 


-143.7(2) 
-145.6(2) 
-147.8(1) 


exp. 




-28.30 


-31.99 


-92.16 


-127.62 
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FIG. 3: (color online) /V max -extrapolated ground-state energies of 4 He 
and I6 as function of the flow parameter a for the NN-only (•), the 
NN+3N-induced (♦), and the NN+3N-full Hamiltonian (a). 



tions for the ground states of 12 C and I6 starting from chiral 
NN+3N interactions. By combining the IT-NCSM with the 
/r-coupled storage scheme for the 3N matrix elements we 
are able to reach model spaces up to A^ max = 12 for the upper 
p-shell at moderate computational cost. Previously, even the 
most extensive NCSM calculations including full 3N interac- 
tions were limited to N max = 8 in this regime 111 611 . As evident 
from the 7Y max -dependence of the ground-state energies, this 
increase in A^ max is vital for obtaining precise extrapolations. 

The general pattern for 12 C and 16 is similar to the light 
p-shell nuclei: The NN-only Hamiltonian exhibits a severe 
ff-dependence indicating sizable induced 3N contributions. 
Their inclusion in the NN+3N-induced Hamiltonian leads to 
ground-state energies that are practically independent of a, 
confirming that induced 4N contributions are irrelevant when 
starting from the NN interaction only. Therefore, the NN +3N- 
induced results can be considered equivalent to a solution for 
the initial NN interaction. The I6 binding energy per nucleon 
of 7.48(4) MeV is in good agreement with a recent coupled- 
cluster A-CCSEXT) result of 7.56 MeV for the 'bare' chiral 
NN interaction lfl7ll . 

In contrast to light nuclei the ground-state energies of 12 C 
and 16 obtained with the NN+3N-full Hamiltonian do show 
a significant a-dependence, as evident from Fig.|2c) and (f). 
The inclusion of the initial chiral 3N interaction does induce 
4N contributions whose omission leads to the ff-dependence. 

A direct comparison of the a-dependence of the extrapo- 
lated ground-state energies for 4 He and 16 is presented in 
Fig. [3] For both nuclei, the NN-only Hamiltonian exhibits 
a sizable variation of the ground-state energies of about 25 
MeV (0.7 MeV) for 16 ( 4 He) in the range from a = 0.04 fm 4 
to 0.16 fm 4 . The inclusion of the induced 3N terms elimi- 
nates this ff-dependence. The inclusion of the initial 3N in- 
teraction again generates an ff-dependence of about 10 MeV 
for le O. Note that the induced 4N (and higher) contributions 
that are needed to compensate the ff-dependence for ie O reach 
about half the size of the total 3N contribution in the SRG- 
transformed Hamiltonian. This is evidence that the hierarchy 
of the many-body forces in chiral EFT may not be preserved 
by the SRG transformation. 
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FIG. 4: (color online) Excitation spectrum for the lowest positive- 
parity states (labelled JnT) in 12 C for the NN-only, the NN+3N- 
induced, and the NN+3N-full Hamiltonian with a = 0.08 fm 4 . 



We have analyzed the role of the different terms of the chi- 
ral 3N interaction at N 2 LO regarding the induced 4N terms. 
By setting the low-energy constants associated with either the 
three-body contact (eg), the one-pion term (co), or the long- 
range two-pion terms (C13 4) to zero and readjusting the re- 
maining Co or ce to reproduce the triton binding energy, we 
find a strong reduction of the a-dependence for the C13 4 = 
interaction. Thus, the long-range two-pion exchange contri- 
bution to the chiral 3N drives the induced 4N contributions. 

Excitation Spectra. In addition to ground-state energies 
the IT-NCSM provides full access to spectra and all spec- 
troscopic observables. As an example we consider the ef- 
fect of 3N terms on the excitation spectrum of 12 C. The a- 
dependence of the excitation energies follows a similar pattern 
as the ground-state energies but is much weaker. The variation 
of the excitation energy of the first excited 2 + state across the 
ff-range discussed before is on the order of a few 100 keV for 
the NN+3N-full Hamiltonian. Thus, meaningful calculations 
of excitation spectra are already possible at the NN+3N-full 
level without accounting for the induced 4N contributions. 

To illustrate the sensitivity of the excitation spectrum on 
the induced and initial 3N interactions we present in Fig.|4]the 
first six excited states of positive parity for fixed a = 0.08 fm 4 . 
Overall, the spectrum agrees well with the calculation pre- 
sented in |01 for Af max = 6 using a different LS -transformed 
chiral 3N interaction. Considering the changes in the spec- 
trum caused by the inclusion of the induced 3N interaction 
(change from NN-only to NN+3N-induced) and the initial 3N 
interaction (change from NN+3N-induced to NN+3N-full), it 
is evident that the induced 3N terms only lead to an over-all 
compression of the spectrum, whereas the initial 3N interac- 
tion affects the individual levels in very different ways, thus 
changing the level ordering. In comparison to the experimen- 
tal spectrum, the agreement for the first 2 + and 4 + states are 
significantly improved by the inclusion of the initial 3N inter- 
action. The first 1 + state, however, is pushed too low by the 
initial 3N interaction, ending up in clear disagreement with 
experiment. The only other state which is not well described 
is the first excited + state, the famous Hoyle state. Though 



the 3N interaction improves the situation, the NCSM basis is 
not well suited for the description of strong a-clustering in 
this state, which gives rise to the very slow convergence (for a 
complementary approach based on chiral EFT see [18)). 



Conclusions. We have presented the first ab initio nuclear 
structure calculations for nuclei in the upper p-shell based on 
SRG-evolved chiral NN+3N interactions. By introducing JT- 
coupled 3N matrix elements and the importance truncation of 
the NCSM model space we are able to surpass all previous 
NCSM studies including 3N interactions regarding particle 
number and model-space size. We find that, unlike for light p- 
shell nuclei, SRG-induced 4N contributions originating from 
the long-range two-pion terms of the chiral 3N interaction are 
sizable in 12 C and 16 0. Based on this finding alternative SRG 
generators can be explored with the aim to suppress induced 
4N terms. Together with the advances presented here, this 
will enable detailed ab initio studies of p- and sd-shell nuclei 
starting from full chiral NN+3N interactions. 
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